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ABSTRACT 



The accuracy of proposed lower confidence limits for system 

availability is analyzed. Random values of the lower 100(1- cO% 

A 

confidence limit L(p<^ ) ^ or s y stem availability are computed for 
a system whose failure density is exponential ( /\ ) and whose repair 
density is exponential The system is modeled as an alternating 



renewal process. The 100 ( 1 )th percentile point of the generated 



distribution of A, 



SL( oC ) com P arec ^ w ith system availability as a 



A 



measure of accuracy for A^ )' 
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I. INTRODUCTION 



The system considered in this thesis is of the type which operate 
for a time, fail, are repaired and returned to operation. The system 
then may be said to have two states and can be modeled as a renewal 
process. 

The objectives of this paper were the following: 

1* To compute exact values of the steady state availability for 
a system with given distributions of times to failure and times to 
repair. 

2. To compare the availability of a system based on the test 
procedures used with the 100(1- th percentile point of the dis- 
tribution of a proposed 100(1 - )% lower confidence limit. This 

will allow an accuracy analysis of a proposed lower confidence limit. 

The method used was to fit a normal distribution to simulated 
data and compute the 100(1 - o( ) th percentile point of the distribution. 

An alternative method was investigated, which fitted a 2 para- 
meter gamma distribution. The procedure proved to be inaccurate. 
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II. STATEMENT OF THE PROBLEM AND ANALYSIS 



A. DEFINITION OF AVAILABILITY 

In this paper the term availability [l] is defined as the probability 
that a system will be n up n or operable when called upon at a random 
time. 

Normally, availability is not specified alone. Both reliability 
and maintainability are usually specified together with availability. 

For durable, continuous operated hardware, reliability can be specified 
by mean-time-to-failure (MTTF). Given the exponential form of time 
to failure distribution, which is applicable to most types of complex 
equipment, MTTF is a constant, the reciprocal of the failure rate. 

Maintainability can be specified by mean-time-to-repair (MTTR) 
and mean-preventive -maintenance -time. Corrective maintenance 
times are often found to be well described by the exponential distribu- 
tion. Given that the time-to-repair is distributed exponential, MTTR 
is a constant, the reciprocal of the repair rate. Accordingly we 
make the following definitions. 

(1) A(t) - Availability at time t, or the probability that 

an item will be operable at a stated instant in 
time . t 7 0 

(2) A^ - Interval availability is the time average of A(t) 

during intervals of length T, and is readily 
obtained from the equation for the average value 
of a function. 
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7" 



/ 



A(t)dt 



*/ 



A(t ) dt 



/ dt 

y o 

c 

Alternately, it may be considered that there is 
some probability distribution h(t) on demand time. 
Then the interval availability is given by 



E [A(t )] 



Cx 

-I 



DO 

A(t) h(t) dt 



In the special case where h(t) is uniform on the 
interval [0, T] 



h(t) = i 



0 ^ t 



o t y t 

X 



thenAp = f J dt ' 

o 

(3) A - the steady- state availability is the limiting interval 

availability as T 

TT 

hf 



A = lim ~ / A(t) dt 

o 



T oO 

It can also be shown that 

A = lim A(t) 

T — > t>o 



If the components do not undergo checkout or repair its interval avail- 
ability approaches zero. For systems with exponentially distributed 
failure times where 



A(t) = 



A 



T 








1 - e 



At 

A— ^ 0 as T~^ taO 
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Thus the steady state or long term uptime ratio of a component which 
does not undergo repair approaches zero. 



If the component is subject to repair, its steady state availability 



is not zero. Assume that both failure and repair times are exponentially 

distributed with means and -y respectively. Given that A(0) = 1 

A A 

i. e. , component available at start of mission, then 

A + ^ c -t(A+ ^ ) 



A(t) = 



A y u -t-A 

yL't*- repair rate 



failure rate 



then A = lim A(t) = — 

t A 

t—> gO 

MTTF 



MTTF + MTTR 



Catron [2] has shown that exact values of A(t) can be and were 

computed for repair distributions other than exponential. The use 

of interval availability as t — ^ i. e. , 

T 

A t = A f A(t) dt A 

o 

to approximate the true availability over the first mission time as is 
currently being done in practice, is a conservative procedure. 

The amount by which the limiting interval availability under- 
estimates the true availability during the first mission time is 
dependent upon the limiting value, i. e. , the lower limiting value, 
the more it underestimates the true availability. 

(4) Ag - System availability. For a series system the 
system availability is defined as: 
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1C 

A c = TT A. 

s •_/ 1 
l>~f 



where A. is the availability of the component of 



type i. i = 1 , 2, 



K 



B. PROBLEM 

The problem considered in this paper is to check the accuracy of 

A 

a proposed lower confidence limit procedure A c _ , » for A c based 

upon component failure time data and repair time data. 

(1) When assuming exponential failure rate ^ 

and exponential repair rate fo . then A. is given by 

(3. 1 1 

A. = 

1 i 

The following ranges of the failure rate ( (?^.) and repair rate 
( '.) has been used in the simulation. 

0. 005 £o(l 0. 01 
3. 0^r I$i4t 10. 0 

which implies that the ranges of MTTF of component i, 
denoted by X . and MTTR of component i, denoted by^C*. 
are as follows. 

loo X i 200 

L lr*£ JL 

1 0 - - 3 

(2) The equation of component availability can be written as 

/&■ _ _ i i 



A. = 

1 /fc+oil 






1 1 + KiMC 

ft* / 

but a. = -r+^r^. = 1 + (<xy*L) 2 - 
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provided the series converges, which it does since 
Since the product q£i is very small, in this case 
0. 0005 £ eC 6 0. 00333. . . 

1 +i \ijZi = 1 ' 

X * 



A. = 

l 



The 



" A s = f ^ 7 r 

Ul i-=\ / L ' 1 y 



Since and JLc' u will have small variance we shall apply the 

central limit theorem, and fit a normal distribution to 



*C 

a s = 1 - Z U 



A 



Then = 1 - 



i.. 

A 

2_oC 

i-l ' 



A A. 



A 



Assuming that A^ 



N 



[ A s- h 



M I H: t I dip, 
2 - ~ Ms 

n£ f To. 






J 



J- 



A A I K '/■/2 ' 1 

implies that A ^ = A £ - Z. -fjf ± 

1 **' X Tu; I 



c< 



f-/ 



C. ALTERNATIVE METHOD 

An alternative method was investigated, which fitted a 2 para- 



meter gamma distribution by method of moments to the distribution of 

A 

, . A 

A 

A, 



A 



A 



<C £ 4C (be 

■ irA ' r ^7* 



6*1 

where^ ^ and^ ^ are the estimates of ^ and ^ ^ respectively. 

Under the assumption that time to failure and time to repair of 
the components were distributed exponentially with failure rate <X • 
and repair rate ^ respectively. The following approach was used 
to make the fit: 
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-L, t . 

0—) I 0 — j I 



Q. = 1-A. = 1 - .. 



<KL 



cic. 'f {3c 

Thus by Taylor's expansion 

- ^ A s = ^L ~ ^(i-Q.) = /L (Q. + ) 

i-l C*i 2 

f ( Qi + $.) « X ; Y; -7$ 

£ = f >■ = £ (ai 4 +u ) 






a. and b. are choosen so that 
i i 

/\ 



c-l 



E( y. - y. 



and 



A 

A a oC i 

Q. = 1 - A. = 1 - 



A 



A. 



A . £ 
Otc * 



where . and . are the estimate of (X . and ^C? ^ respectively. 
The 2 parameter gamma, r<r, V ) was fitted by method of 

A\ 

moments to the distribution of • 

E( / 5 ) = — 

v"- < y s > = -71 

A 1 

e ( A ) . y" 



Therefore 



A 

2 - y« > 



E { 



Then 



v .u- ( £) 

A 

r = 



= r 




far (Y*) 
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S' 



by replacing all oC . and ^3 ^ by &C ^ and ^ 
Then 2 Y Y 6 is distributed 

C 



sir 



s 

A s = e 



r_ 

r 



Y 



1 . e. 



From this assumption a one sided lower 100(1 - )% C. L. for A c 

i. 

was computed. 

P(2YY 6 >% <)2r ) = 1 - o< 

P( e~y Z Ctf> ^ li r ^- \\ = | — ^ 

p(a q > f)sL Lot-) = i - <*: 



r i 

where L 2 r j = smallest integer greater than or equal to 2 r 
The procedure showed to be inaccurate due to small values of p 



k: 



A- /V 



D. PROCEDURE 

The procedure used to fit a normal distribution to 1- jLc<~ 

f 

as described in B is based on computer simulation. 

The computations of the estimates ofO^^ and s (1 = 1> 2, ' • # k) 

are based on the assumption that we are given a series system with 
K types of components. 

V.j - time to failure of the j-th component of type i. 

The time to failure is assumed exponential with 
failure rate (X .« 

T^. - time to repair of the j-th component of type i. The time 

to repair is assumed exponential with repair rate 

V.. and T.. are independent. 
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(1) A random sample of size 1NL components of type i are tested 
until failure or a specified planned test time (PTT). 

is the PTT of the j-th component of type i 



X.. = J 
ij ' 



i if the j-th component of type i fail, i. e. , 

life time less than T .. 

oij 



O 



if the j-th component of type i lives beyond 

T ... 
oij 



>L 

/ X c. i* = # of failures of type i components. 

Jr I J 

The components were tested against the (PTT) T q .. and the observed 

operating times V .. set equal to V. . or T .. whichever is the smaller. 
r oij ij oij 

(2) A random sample of size n. of failed components of type i 

are repaired with repair times T.^, T.^, . . . , T.^ . 



(3) Using the test data described in C 1 and 2 the estimates of 

(X • andAU* were computed. 

1 / 1 * 

a At 



A, 



A 



A/*; . 

2 14 li 

& — 

• = t ' <■* 



2 hJi +■ I 



j; 



=i 



tli 



(4) With the data in C-3 1000 values of A g ^ ^ ^ were computed. 



A A 

A SL(« ) = A s 



/ ^ 
ft 



X \t_ 

^ Tel: I A 



j-l J 
A KL ^ ^ 

where A g = 1 - X oC C LX_C 

i-t / 



and 
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z 

c* 



is the 100(1- ) th percentile point of the N(0, 1) distribution. 



E. ANALYSIS 



A 



The generated values of Ag ^ ^ ^ was ordered from lowest to 

highest and the 1000(1- OC)% one was compared to Ag. 

A 

The Ag ^ was a-l so compared to Ag and the probability of 

success was computed by the formula 

A 



P (A S L( <X ) 



- A e ) 



based on the 1000 values of A, 



S L( G< )' 



The results are shown in Appendix A for the different combinations 

of , N., n. and T ... To avoid fluctuations due to the 

i 1 ° i i i oij 

random number generator, the values of Ag ^ ^ ^ at different^ 1 s 
was computed by only changing Z^ 
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III. ACCURACY RESULTS 



By changing the values of the input parameters the accuracy of the 
procedure was checked. The following conclusions were made based 
on the simulation: 

a. Use of the steady state availability of the system (Ag) to ap- 

A 

proximate the estimated availability Ag ^ ^ j over the mission time 
is a conservative procedure. 

b. The amount by which the steady state availability differ from 
the estimated 100(1- )th percentile point is very small. 



Example 

Case 3 - Appendix A. 

With the given parameters the systems steady state availability is 
computed, giving 



kL 

A„= T 






b -I /3c + 



= .9925 



The computed 95th percentile point of the distribution of Ag 
is equal to . 9931 

. 0006 



Thus, 



A 

A S “ A S(950) 



9925 - . 9931 



Similarly the 90th and 80th percentile points of the distribution of 
A A 

A S L( 1) anC * A S L( 2) are ' ^930 and ' 9929 respectively. By increasing 
the number of items which are life tested from 50 to 100 and keeping 
the other parameters constant, we get the result given by case 5. 

That is, 

A„ = 



A 

A 



S L(0. 05) 



9925 (as before) 
= . 9928 
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This implies that 



A 

A S " A S( 950) 
A 



. 9925 - . 9928 = 



0003 



A 

In this case 90 O) an< ^ A S( 800) = * ^926 w ^ich implies that 



A 

A„ - A 



SL(oC ) 



= . 0001 



for 0C = 0. 1 and = 0. 2 



The true levels of confidence associated with A . . as a lower 

o -L/( Ox ) 

100(1- (X ) confidence limit for A are given in columns 10 through 12 

S A 

for = • 05, . 10, . 20 respectively. That is for Case 3, Ag ^ 

is an 90, 8% lower confidence limit for Ag rather than 95% lower 

A 

confidence limit. This is a measure of the inaccuracy of A 0 T . 

o L(, Uo ) 

as a lower 95% confidence limit for A Q for the given parameter values 

A 

in Case 3. Likewise Ag ^ is really a 84. 5% lower confidence 
A 

limit and Ag ^ i s really a 76. 9% lower confidence limit. 

This variation between the true level of confidence and the proposed 

level of confidence (e. g. , 90. 8% vs. 95%) are due to the very small 

A 

variance of A T , > * for the cases considered. It is felt that for 

b Lt{ Ox ) 



these cases the quantities 



A 

*S " A s(l- <X ) 

are better measures for the accuracy of this procedure. 
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Result of the Simulation 
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APPENDIX B 



LEVEL in 



’A I M 



DATE = 710A6 



09/A 1/50 



OTHERS T CM_ ALf A 1 1 5 ) , b LTA( 15) , V ( X5 . 1 00 ) * X ( 1 5 , 1 0 C ) * T ( lb,60J,- 
l AM AT ( l 5 ) , YHAT ( 1 5> ) , U( 15 ) , M I U 5 ) , AS 1 ( 1 000 ) » AS2 ( 1000) , ' B 3 < 1000) 

,-lFAL PTT ( 1 5 , IOC ) /I 500*2 .0/ 

f-0\= 1 ( 00 



K - 1 5 
Z 05= 1 



6A5 



100 



1000 



5(0 
6( C 



zero 



Z 10= 1.20? 

7 20= C. M A 2 
I X=o 3A5 13 .. 

T C v MAT ( 2 F 1 o 
) COO 1=1 



5,2110) 

, < 



RCjc (3,100) «LFA( I ) ,3ETA< I ) ,NU ) .NI< I ) 



CCMT I MJE 
X 1=0 . 0 
X 2=0.0 
X T = 0 . 0 

NRlTEt 3,500) ----- 

F r t- ! M A T ( / / / . 3 T X , 1 ASHA T 0 5 * 
f 0 KM AT ( 30X , 3F10.A) 

A S = 1 .0 

OC 7( CO 1 = l , K 

A$ = AS*4F-TA( 1 )/ ULF A ( I ) 4-bETA( I ) ) 

nc coco ;•'=), on 
no 2000 l=l,K 

i\\ = N ( 1 ) 

CO K 00 J = 1 , r N 

X ( I , J ) =0.0 

CALL 0 AMOU ( I X , IV ,YFL ) 

1 X= l Y 

, 0 / A L F A ( 1 ) *ALO 

,IT.PTT( 1 , J) ) 

, b T . 0 T 7 ( I , J ) ) 



3 X , * A S H A T 1 C ’ » 3 X , 'ASHA 120' ,///) 



V ( I , J ) =- 1 
1 F ( V I I , J ) 

I F ( V ( I , J ) 

JCOO CCUTIMUE — 

?oo(3 r.CMriMjr 

no 3 x? a i = i , x 

h .\ T = M ( I ) 

f:° , A A A J = 1,CMI 

C ALL « At.UU t I X , I V , V r L ) 

T( T ! J)=— l . C V L> 0 T A ( I ) *AL n G( VOL ) 

_A A A A CC.NTINJF - • — 

13’? CONTINUE 
V A=0 .0 
A YF A T= o . 0 
!)G ACCO I =1 , •< 

SU V T = 0.0 

SU v XX=0.0 

SLYV=e.C 

_ - .\N = M I ) . ... — 

no 5CCU J=1,NN 

S L M X X = SU M X X +• X ( I , J ) 

S U v V = S u , J V 4- V ( I , J ) 

SU YT = SUi* T + P T T ( I , J) 

3)00 CCM INUF 
SU V TT=C . 0 
N M = M ( I ) 

DO 5 5 55 J = 1 ,NM — 

SL' 1 T T= S >J M TT fT ( I , J) 

5 5 r 5 CC NT INUF 

X N = M ( I ) 

X M = M I { I ) 

A FAT ( I ) = SU' 1 X X / S I'YV * ( 2 
YHAT ( T )=SU"T7/XNI 
V A = V A + CHAT ( I ) * Y F AT ( I ) ** MSUOT 

A YHA T = A YHAT 4- AHA I I I ) ’c YHA T ( I.) 

A 000 CONTINUE 

A SH A T = 1 .O-AYFAT 
A S 1 ( F ) = A S H A T- S OR T ( V A ) * 7 05 
AS2 ( Y ) =ASHAT-SCRT( VA ) U 'Z 10 
ASM F ) =ASHA T-SOkT{ VA ) * 7 2 n 
I F ( A S .OF . AS 1 ( M ) ) X 1 = X 1 + 1 . 0 
I F ( A S . OF . AS 2 ( '•■ ) ) X 2= X 2 4- l . 0 
iriAS.GI .AS3(.'L)) X 1 =-X3 *-l . 0 



MY FI. ) 

X ( I , J ) = 1 . C 
V(I,J)=PFT(1,J) 



0 *XN ) / (2 .0*XN+ l .0 ) 
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LFVEL 18 MAIN DATE = 71046 09/41/50 

8 J C 0 CONTINUE 

DC 8800 1=1,599 
J J = L + 1 

DO 8°fO J=JJ,1000 
I E ( AS1 ( L) . L r .ASl (J) ) GO TO 800 
T F P P 1 =AS 1 (L ) 

ASl ( L) =AS1 ( J ) 

ASl ( J) =7 E'-iPl 

SCO J F ( AS2 ( L ) . LT . AS2( JJO—GC— T-CL_aSQ 

TFPP2=AS2 ( L ) 

A S 2 ( L ) = A S 2 ( J ) 

A S2 ( J 1 = T E M P 2 

RPO 1 F ( A S 3 ( L ) . L T . A S 3 ( J I ) GO TO 3800 
TEMP 2= AS 3 ( L ) 

A S3 ( L ) =AS3 ( J ) 

A S 3 ( J ) = T F V P 3 

3800 CONTINUE .. .... 

WRITE (6,600) (ASK I ) , AS2( I ) ,AS3( I ) , 1 = 1, 1000) 

WRITE (6, POO) 

200 FORMAT <// / ,32X , ’ A L r A ( I ) ’ , 3 X , 1 BETA( 1 ) * , 7 X , * N ( I ) 1 ,5X, • N N ( I ) * , ///) 

DC J = l , K 

W R I T 0 ( 6 , 3 00 ) ALF A ( I ) , 8 E T A ( I ) , N (I) , M ( I) 

9 9 C 9 CONTINUE - ■ 

300 FORMAT (30X,2FIO. 5, 2110) 

P R 9 3 1 = X 1 /X'”l 
PR0ii2= < 2 / X 
p r ■>, - < 7 / v •• \. 

.Oil (• ! O, V ) AS , P .031 , PRO 82 , » J P 0 0 3 , M\ 

4 00 F 0 R M A r ( / / / , 1 0 a , 1 A S = ’ , P 7 . * , 4 X » * PkiJUI = ' , F 7 . 4 , 4 X , • PR C 3 2 = • , F 7 . 4 , A < , 

1 • PRO 8 3 = ' , F 7 . 4 , A X , * ’N = • , I 4 ) 

w 0 I T F ( 6 , 00 ) AS 1 ( : bO ) , AS? ( 5 00) , A S 3 ( 3 • 1 0 ) 

700 FOR" AT (// / , 1 OX , ‘ AS 1 ( ^'5 0 ) = ' ♦ T 7 . 4 , 4 X , ' A S 2 ( 900 ) = ' , F 7 . 4 , 4 X , ' AS 3 ( 300 ) = ' 

2 , F 7.4) 

SI CP 
FMD 



APPENDIX C 



A 

Derivation of Mean and Variance of A r 



K /I M [ 

A s -r Hi -it / o.,>. T i-i- ^ 

b i~/ ^ *- l '~ ( 



/ 
^3^ 



i** 



but -T — jU~1 = MTTR 



The MTTR should always be less than one mission unit, which implies 
that 

< ^Ai < ^ i 



thus 



1 



A. = ^ . 

1 1 -f a 






= i -0^1+ (ocj ^ ) 2 - (^; ) 3 + 



provided the series converges which it does since i— 

Since the product C\i.jLci^ s very small 



A i 1+ 1 



OCCft^ 



Thus 



A s = II 



tc 

If 

c=l 



ic 1 



1“^'^ 

Since £*C ^ and ^ ^ will have small variance we shall apply the central 



limit theorem, and fit a normal distribution to 

AL 

^ 0^ c LsC-L* 

0~t 



A s = 1 



Thus 



A J6 A A 

A g = 1 - ^ Ac 



A A 

when^/ ^ is the estimated MTTR and C>d ^ is the estimated failure 

rate. Based on the assumption that time to failure and time to 



19 



repair is distributed exponentially we can find the mean and variance 

A 

of Ag. 



A 

E(A ( 



, ) = i -Zfc(*0 kdh — I iLct- iu- i — & 
> * / /*%»/' 



A 



By using the equations given by OD 29304 [3], for^ . and Vctr (oi .) 

A 

we can find the variance for Ag 



A 



. _ # of failures 
# of test time 



2 Ni 



A 



2 N. 4 1 
i 



E( X i ) = 0 < 

A 

VOX- ( . ) = 

l 



oC c 



Ml 

where 2 £< 






j 



A 



= sum of all planned test times on component 
i (in mission units). 

ru 

2. / c j 



//Lei - ft repair times 



n. 

i 



J 



I-/ 



n L 



z 



A 



E (^x i ) 

VOS-(llC) -- /±± 
/ ru 



U /V 



Thus VOA- (Ag) = V (5U~" (1 - t <* i ) 



Vur/ 






But 



- t(£c % )£yi) -oil/u 

E(o< •) = V*W- (<£,;)■+ ^ 



( 1 ) 



2 

c 



cU ' , a 
+ C< 0 

■Z ti-. 

J J 
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E (fat) = l^r (fa ) •+ fat 

V“* ( 



Therefore 



fafafa^i^-. +*<. Xfa'i^r) - °V 

i 



2 ft;*/ \ 



2 



'2. 



__ oCt JU j* (li-f-l , .2 hi '^1 ^- 



+ <* 



T id n ’ 



A 



Hw 



which implies that 
/v 



vo-r 



(X .<_ 'j 



* °( ^ Hi +1 



Z Aj n. 



Thus by ( 1 ) 

M A \ . -r n £ ■+/ 

var (/f-s) - f iv; -— 

^ ’Z ii.- n - 

j*» J 

A A 

Since X . and is approximately equal to X, ^ and^/x. . respectively, 



and 



n. + 1 1 

i 1 a s n . — oo 



n. 

l 



we can write 

A 

L s 



Jfa /\ /v 

E(A q ) = / — "2. Oi C jU-o 
c— ^ / 

a *r J: 

VOS' (A„) = > " 

s ^7? 



, = /? 



c-/ 



Zt L 



-/ 



J 
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which implies that 



A 

A g ^ N(A S , 



Var (Ag)) = N(Ag, Var (Ag)) 



= N 







2 . 



A A 2 

r.uj 



Thus by the above assumptions we can write 

p ( A s' A s - ^) ~ /- 

I /t Vr 

7 ( £ 4 - /?> £ ij'MC)) = /- o< 

T L~ ft* ~ I ~ 

? t /?-> > /? s - ^ \CTr) = 1 - ~ 



A 

Thus A 



SL('X ) 



A 

= A S ■ Z <A 



{Yar (/} s) 



= A s - Z c< 



<*i /Y2~ 



i -~l 



K' 

1L tdj 

i 



1 

2 - 
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